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We investigate heat and charge transport through a diffusive SIF1F2N tunnel junction, where N (S) is a normal 
(superconducting) electrode, I is an insulator layer and Fi 2 are two ferromagnets with arbitrary direction of 
magnetization. The flow of an electric current in such structures at subgap bias is accompanied by a heat 
transfer from the normal metal into the superconductor, which enables refrigeration of electrons in the normal 
metal. We demonstrate that the refrigeration efficiency depends on the strength of the ferromagnetic exchange 
field h and the angle a between the magnetizations of the two F layers. As expected, for values of h much larger 
than the superconducting order parameter A, the proximity effect is suppressed and the efficiency of refrigeration 
increases with respect to a NIS junction. However, forh^A the cooling power (i.e. the heat flow out of the 
normal metal reservoir) has a non-monotonic behavior as a function of h showing a minimum aih^ A. We also 
determine the dependence of the cooling power on the lengths of the ferromagnetic layers, the bias voltage, the 
temperature, the transmission of the tunneling barrier and the magnetization misalignment angle a. 

PACS numbers: 74.45.+C, 74.50.+r, 74.25.fc, 75.30.Et 



I. INTRODUCTION 

The presence of the superconducting energy gap leads to 
a selective tunneling of high-energy quasiparticles out of the 
normal metal in a normal metal - insulator - superconductor 
(NIS) tunnel junction.^ ^ This phenomenon generates a heat 
current from the normal metal to the superconductor (also re- 
ferred to as "cooling power"). The heat transfer through NIS 
junctions can be used for the realization of microcoolers.'^'^ 
Present state-of-the-art experiments allow the reduction of the 
electron temperature in a normal metal lead from 300 to about 
100 mK, offering perspectives for on-chip cooling of nano- 
sized systems, such as high- sensitive detectors and quantum 
devices."^ 

The cooling power of tunnel junctions depends on several 
parameters, some of them controllable. For example the cool- 
ing power can be optimized by controlling the voltage across 
the junction. A maximized cooling effect is reached at a volt- 
age bias just below the superconducting energy gap A. Larger 
values of voltage, eV > A, lead to a larger charge current / 
through the junction and hence to larger values of the Joule 
heating power, i.e. to a negative cooling power. A limitation 
of the performance of a NIS microcooler arises also from the 
fact that nonequilibrium quasiparticles injected into the super- 
conducting electrode accumulate near the tunnel interface.^ ^ 
As a consequence hot quasiparticles may tunnel back into the 
normal metal, leading to a reduction of the cooling effect.l^^ 
In order to overcome this problem a so called quasiparticle 
trap,^^ made of an additional normal metal layer has been at- 
tached to the superconducting electrode, removing hot quasi- 
particles from the superconductor. Recently, it was also shown 
that a small magnetic field enhances relaxation processes in a 
superconductor and leads to significant improvement of the 
cooling power in NIS junctions. ^^ Improved cooling perfor- 



mance can be also achieved by proper tuning of the tunneling 
resistances of the individual NIS tunnel junctions in a double 
junction SINIS cooling device.^ 

Another important limitation for NIS microcoolers arises 
from the intrinsic multi-particle nature of current transport in 
NIS junctions which is governed not only by single-particle 
tunneling but also by two-particle processes due to the An- 
dreev reflection. ^^ While the single-particle current and the 
associated heat current are due to quasiparticles with energies 
£ > A, at low temperatures or high junction transparencies the 
charge transport in NIS junctions is dominated by the Andreev 
reflection, i.e. by subgap processes. The Andreev current Ia 
does not transfer heat through the NS interface but rather gen- 
erates the Joule heating IaV. At low enough temperatures this 
heating exceeds the single-particle cooling. ^^^^ The interplay 
between the single-particle tunneling and Andreev reflection 
sets a limiting temperature for the refrigeration Tmin- 

One way to decrease T^nm is to decrease the NIS junction 
transparency. However, large values of the contact resistance 
hinder carrier transfer and lead to a severe limitation in the 
achievable cooling powers. In order to increase the junction 
transparency and at the same time to reduce the Andreev cur- 
rent, it was suggested to use materials where the proximity 
effect is suppressed, such as ferromagnets, ferromagnetic in- 
sulators, and half-metals. In particular Giazotto et at. stud- 
ied theoretically a ballistic normal metal - ferromagnet - su- 
perconductor structure within a phenomenological model and 
predicted an enhancement of the cooling efficiency compared 
to NIS junctions. ^^ The reason for that increase lies in the sup- 
pression of the Andreev reflection due to the band structure 
of the ferromagnetic metals. The electron involved in An- 
dreev reflection and its time-reversed counterpart (hole) must 
belong to opposite spin bands; thus, suppression of the An- 
dreev current occurs in a FS junction and its intensity depends 



on the degree of the electron polarization at the Fermi level 
which is proportional to the exchange field of the F layerP^^*^ 
The enhancement of the cooling efficiency by the magnetic- 
field-driven tunable suppression of the Andreev reflection in 
superconductor/two-dimensional electron gas nanostructures 
was also studied in Ref. |23] Note that theoretical studies 
of electron cooling in S F prox imity systems were performed 
only in the ballistic caseP^EIl while real metallic systems are 
in the diffusive limit. Moreover, ferromagnets show in general 
a multi-domain structure that was not considered in previous 
articles. 

In this work we present a quantitative analysis of the ther- 
moelectric transport in NIS microcoolers with a diffusive fer- 
romagnetic interlayer consisting of two magnetic domains 
with arbitrary direction of magnetization (so called "super- 
conducting triplet spin-valve'^^). Based on the quasiclassical 
Keldysh Green functions formalism we compute the electric 
and heat currents through the junction. We show that the en- 
hancement of the cooling power with respect to the NIS case, 
as proposed in Ref. 18 , only works if the exchange field of the 
ferromagnetic interlayer h is much larger than the supercon- 
ducting order parameter A. However, the cooling power shows 
a minimum value for h^ A. We also study the dependence of 
the cooling power on the angle a between the magnetizations 
of the ferromagnetic domains. In the case of weak ferromag- 
nets, i.e. for /z < A, the antiparallel configuration a = n leads 
to higher values of the cooling power and smaller values of 
Tmin- By increasing h this behavior is reversed and the mono- 
domain configuration (a = 0) is more favorable for the refrig- 
eration. For large values of the exchange field, h :^ A, the 
cooling power is almost independent on a. 

The paper is organized as follows. In the next section, 
we formulate the theoretical model and basic equations. In 
particular, we obtain the expressions for the electric and the 
heat current and identify the contributions corresponding to 
single-particle and Andreev tunneling events. In section III 
we present and discuss the main results of our work. We fi- 
nally conclude with a summary of them in section IV. 



II. MODEL AND BASIC EQUATIONS 

We consider the SIF1F2N junction depicted in Fig. [T] A 
ferromagnetic bilayer F1F2 of length /12 = /i + /2 smaller than 
the inelastic relaxation lengtlP^ is connected to a supercon- 
ductor (S) and a normal (N) reservoirs along the x direction. 
The F1F2 bilayer can either model a two domain ferromagnet 
or an artificial hybrid magnetic structure. We consider the dif- 
fusive limit, i.e the elastic scattering length i <C min((^/^,(^), 
where ^h = ^/WJ^ i^ the characteristic penetration length of 
the superconducting condensate in to the f erromagnet, h is the 
value of the exchange field, ^ = ^/WJ2A is the superconduct- 
ing coherence length and & is the diffusion coefficient (we set 
h = ks = I and for simplicity we assume the same & in the 
whole structure). 

We also assume that the F1F2 and F2N interfaces are trans- 
parent, while the SFi is a tunnel interface. Thus, the two ferro- 
magnetic layers are kept at the same potential as the voltage- 
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FIG. 1. The SIF1F2N junction. The interface at x = corresponds 
to the insulating barrier (thick black line). Interfaces ai x = li and 
X = li2 are fully transparent, a is the relative angle between the 
magnetization directions of Fi and F2. 



biased normal reservoir. The magnetization of the Fi layer is 
along the z direction, while the magnetization of the F2 layer 
forms an angle a with the one of the layer Fi . Both mag- 
netization vectors lie in the yz plane. Correspondingly, the 
exchange field vector in the Fi is given by h = (0,0,/z), and 
in the F2 layer by h = (0,/zsina,/zcosa), where the angle a 
takes values from (parallel configuration) to n (antiparallel 
configuration). 

In order to describe the heat and electric currents through 
the structu re we introduce the quasiclassical matrix Green 
function gW^ 



G^ G^ 
G^ 



(1) 



The latter is a matrix in the Keldysh x Nambu x spin space. 
The R, A and K indices stand for the retarded, advanced and 
Keldysh components (we use the symbols " for 8 x 8 and T for 
4x4 matrices). By neglecting non-equilibrium effects, the 
Keldysh component is related to the retarded and advanced 
ones by 



&h- 



fiG^, 



"^ = H'"'^ 



n - 
eV 



- n^ + T^H- 



±tanh- 



E-eV 

2Tn 



(2a) 
(2b) 



where n± and T^ are correspondingly the equilibrium quasi- 
particle distribution functions and the temperature in the nor- 
mal reservoir and T^ is the Pauli matrix in Nambu space. 
The retarded and advanced components are related via G"^ = 

The matrix Eq. ([l|) obeys the Usadel equation,^^ which in 
the notations of Ref. 1 3 Oi reads 

i^dj = [t, {E - hcj) , G] , / = Gd,G, G^ = 1 , (3) 

where cr = (dx, Oy, cr^) are the Pauli matrices in spin space. In 
the Fi region ha = ha^ and the Usadel equation Eq. ^ has 



the form 



i9dJ=[x,{E-a,h),G\, g2 = 1. 



(4) 



In the F2 region her = ha^&x^{—i(yxa). It is convenient to 
introduce Green's functions rotated in spin-spacep3 



G = U^GU, U = exp (ia^a/l) . 



(5) 



The rotated function G is then determined by Eq. (|4|). 

The Usadel equation Eq. ^ should be complemented by 
boundary conditions at the interfaces. As mentioned above, 
we assume that the F1F2 and F2N interfaces are transparent 
and therefore the boundary conditions at x = /i , /12 read 



lx=/i-0 ■ 



dM 



x=U -0 ■ 



d,G 



x=/i+0' 



lx=/i; 



-0 ~ ^^' 



(6) 

(7) 
(8) 



At X = 0, the SFi interface is a tunnel barrier, and we may 
use the Kupriyanov-Lukichev boundary conditions,^ 



4=0=(^/^)[^^'^]x=0' 



(9) 



where Gs is the Green function of a bulk BCS superconductor 
defined as 



Gs = t^U + TjV, 



(w,v) = {E, /A)/e, £ = ^J{E^i^Y-^^, (10b) 

where r] describes inelastic scattering rate and W <C 1 is the 
diffusive transparency parameter,^^ W = ^ /'^8nR- Notice that 
W is temperature-dependent, since the coherence length ^ is 
proportional to A~^'^{T). In section III we set rf ::^ 10~^Ao in 
our calculations, where Aq is the superconducting gap at T = 
0. In the following we omit rf in our analytical expressions 
for simplicity. 

Because of the low transparency of the SFi barrier, the 
proximity effect is weak and the retarded Green function can 
be linearized (we omit the superscript R), 



-^x/, 



(11) 



where / is the 2 x 2 anomalous Green function in the spin 
space (1/1 <C 1) that obeys the linearized Usadel equation, 

/^4/ = 2^/-{hcT,/}, (12) 

where {•, •} stands for the anticommutator. 
The general solution of Eq. ([12]) has the form 

/(X) =/(x)+/,(x)CT,+/,(x)CT„ (13) 

where / is the singlet component and f^, fy are the triplet 
components with respectively zero and ±1 projections on the 
spin quantization axis (we choose the z-axis). The fy term is 
usually known as the long range triplet component because it 
describes Cooper pairs with parallel spins which survive the 



strong exchange splitting and can diffuse into the ferromagnet 
over larger distances compared to the singlet component. ^^^^' 
Indeed, by substituting Eq. JT3|) into Eq. ^V2\ and using the 
boundary conditions Eqs. ( [6p] ) one can compute the com- 
ponents /(x), fy{x) and fz{x) , and easily show that while 
f{x) and fz{x) decay into the ferromagnet over the magnetic 
length ^J Q> jlh the long-ranged component fy {x) decays over 
the length given by ^/WJ2E. 

The charge and energy currents, / and Q respectively, can 
be obtained fro m^^ i ^^'^^'^Q' 



/ = 



^-^Tl^dE, Q=^4rEI^dE, 

e Jo e^ h 

(l/8)TrT,y^, /+ = (l/8)TrTo/^, 



J^ = (Gd^G)" = G^d^G" + Cfd,&. 



(14a) 

(14b) 
(14c) 



where /^ is the Keldysh component of the matrix current de- 
fined in Eq. (jSland Tq is the unitary matrix in Nambu space. 
From Eqs. ^9]), ( |lQa| ) and ([13]) we determine both currents, 
/ and Q at the SFi interface 

1 r 
/=— / ^_(Rew + RevRe/o) (i£, (15a) 

eR j{) 

1 r 
2=-^—/ £'(^+ — ^) (Rew — Imvlm/o) J^", (15b) 
e^R J^ 

where n = tanh(£^/275') and Ts are the equilibrium quasipar- 
ticle distribution function and temperature in the supercon- 
(10a) ducting reservoir respectively. The function /o = f\^^Q is the 



singlet component of / at x = 0. Notice that only the singlet 
component of / enters the equations for the electric and en- 
ergy currents. There is however an indirect dependence of the 
currents on the triplet component since the amplitude of /o in 
turn depends on the amplitudes of the triplet fy and f^. 



We now examine Eqs. ( 15 ) and discuss separately two main 
contributions for the currents / and Q, which originate from 
the single-particle (E > A) and Andreev (0 < £^ < A) pro- 
cesses. Let us focus first on the single particle contributions 
/i and 2i. For energies larger than the superconducting gap 



(E > A) only the terms proportional to Rew in Eqs. ( 15 ) are 
non-zero. From Eq. ( |10b| ) and Eq. ( |15a| ) we obtain the single- 
particle contribution to the electric current 



h = 



^ rNs{E)n-{E)dE, (16) 

eR Ja 



whQYQNs{E) = \E\e{\E\ -A)/VE^-A^ is the BCS normal- 
ized density of states (DOS) and 0(x) is the Heaviside step 
function. 

Rewriting in Eq. (16) the n-{E) function in terms of 
the Fermi function in the N reservoir, nF{E) = [1 + 
exp(£^/7/v)]~^ we arrive at the well known expression for the 
tunneling current,^ 

h = ^J^^Ns{E)[nF{E-eV)-nF{E)]dE. (17) 

Note that within the linear approach the normalized DOS in 
the F layer is equal to unity and therefore the single particle 
electric current is independent of /q. 



The single-particle contribution to the energy current can 



be obtained from Eq. ( |15b| ), 
Qx = ^J^E{n^-n) [Ns{E)-M^{E)lmfo] dE, (18) 



whereM+(^)=A0(|^|-A)/V^2_^2^ 

For energies £" < A the electric charge is transferred by 
means of the Andreev reflection. The subgap current or An- 
dreev current can be obtained from Eq. ( |15a| ), 

Ia = ^ [ n-{E)M-{E)RcfodE. (19) 

whereM-(^) = A0(A- \E\)/VA^-E^^ 

According to Eq. ( |15b| ) the contribution of the Andreev pro- 
cesses to the energy current vanishes, Qa = 0. 

We are interested here in the cooling power P, i.e. in the 
heat current flowing out of the normal metal reservoir. One 
can express the cooling power in terms of the contributions 
introduced previously, ^^ 



-Q-IV = Pi^Pa. 



where 



Pi = -Qi-hV, Pa 



-IaV. 



(20) 



(21) 



From Eqs. ( 17|18|19| ) it is clear that for equal temperature of 
the electrodes (Tf^ = Ts) and no bias voltage the cooling power 
vanishes. For a finite voltage V, on one hand the heat is taken 
from the N reservoir and is released in the superconductor. 
On the other hand there is a global heat production in both 
electrodes due to the Joule heating. In particular the Andreev 
current I a contributes to the Joule heating IaV, which is fully 
released in the normal metal electrode and leads to a reduction 
of the cooling power. 

In the next section we calculate the cooling power of the 
SIF1F2N junction as a function of the different parameters by 
solving Eqs. ( |17|18|19|2Q| ). 



III. RESULTS AND DISCUSSION 



According to Eqs. (15 ), in the linear case both the electric 
and heat currents are determined by the singlet component /o 
of the anomalous Green function, E q. (p3] ), evaluated at the 
SFi interface (x = 0). Solving Eq. (|12|) in the Fi layer we 



obtain for the components of Eq. (13 ), 



2W 
f± {^) = ^± cosh(^±x) + - — {ua± — v) sinh(^±x) , (22a) 
k± 



2W 

fy{x) = ayCO^\\{kyX) + " — uUy ^ivih{kyX) ^ 



(22b) 



where f±= f±fz, at are the boundary values of /J- at x = (/ 
stands for +, — , j) and 



*, = ,/?(£^, *. 
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FIG. 2. (Color online) Cooling power versus exchange field for dif- 
ferent orientations of the exchange field vector in the second ferro- 
magnetic layer F2: a = (black solid line), a = Tijl (blue dashed 
line) and a = n (red dash-dotted line), calculated at optimum bias; 
W = lx 10-^ T = 0.25Ao, li=^ and l2 = 6^. We have defined 



yoptj 



'^o/Ag. 



In the F2 layer the general solution has the form, 

fi{x) =bismh[ki{x-h2)], 



(24) 



where fi are the components of the rotated Green function, 
Eq.Q. 

Using the boundary conditions at the F1F2 interface, Eqs. 
([6][7]) one obtain a set of six linear equations for the six coeffi- 
cients at and bt, that can be solved straightforwardly. 

In particular we are interested in the value of the singlet 
component of the anomalous Green function at x = which 
is given by /o = (a+ H-a_)/2. The analytical expression 
has lengthy awkward form and we do not present it here. 
Once we obtain /o we compute the charge and energy cur- 
rents from Eqs. (16), ([18]) and ([19]). Finally, using Eq. (20) 
we determine the cooling power. In what follows we assume 
that the temperatures of the S and N reservoirs to be equal, 
Ts = Tn = T, and neglect nonequilibrium effects in the ferro- 
magnetic interlayer.Q3 

The bias voltage between the S and N reservoirs is an eas- 
ily adjustable experimental parameter, so all our curves except 
those presented in Fig.|5]are calculated for optimal value of the 
voltage bias Vopt, at which the cooling power reaches its max- 
imum for given values of the other parameters. In what fol- 
lows, we assume the quantity W to be taken at T = 0, allowing 
for its temperature dependence in Eqs. ( [22] ) by means of cor- 
responding temperature-dependent factors. In the subsequent 
analysis the cooling power P is given in units of Al/e^Ro, 
y^^) where Aq is the value of A at zero temperature and Rq is the 
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FIG. 3. (Color online) The Andreev current as a function of (a) the 
exchange field for I2 = 6^ and as a function of (b) the F2 length for 
h = 0.7A(r). Different magnetic configurations are chosen: a = 
(solid black line), a = n/2 (dashed blue line), a = n (dash-dotted 
red line). The Andreev current is calculated at optimal bias; W = 
7 X 10-^ T = 0.25Ao, h=^. Wehavedefined/A =/A(Vi,p?)^/?o/Ao. 



junction resistance at a fixed value W = 10 ^ of the tunneling 
parameter. 

We first study the dependence of the cooling power on the 
strength of the exchange field h. This dependence is shown 
in Fig. [2] for three different angles a = 0, 7r/2, n between the 
magnetizations of Fi and F2 layers at the optimum value of 
bias voltage. We have chosen the values of the temperature 
and tunneling parameter W such that the Andreev current role 
in the cooling processes is essential (see Fig. [6]).^^ The thick- 
ness of the F layers is chosen to be /i = (^ and l2 = 6^. 

Depending on the value of h/^h^ where ^h = y^Wj2h is the 
characteristic penetration length of the superconducting con- 
densate into Fi, one identifies different behaviors. If h ^ t,h, 
i.e. for large values of h/A{T) the amplitude of the super- 
conducting condensate in F2 can be neglected, as well as the 
dependence of the /o function on the angle a. Thus, in the 
limit h/A{T) ::^ 1, the value of the cooling power does not de- 
pend on a. Moreover, this asymptotic value is larger than in 
the nonmagnetic case (h = 0). This is a consequence of the 
strong suppression of the singlet correlations in Fi due to the 
exchange field and hence of the Joule heating associated to 
the Andreev current [see Eq. ( [2Q| ] . Note that for the value of 
temperature used in our figures A{T) ^ Aq . 

In the opposite limit, /i/(^/j <C 1, the characteristic penetra- 
tion length depends weakly on h, and therefore the cooling 
power is also a-independent. However, by increasing h the 
cooling power first decreases and reaches a minimum. This 
unexpected behavior is qualitatively similar for all magnetic 
configurations and is a consequence of the Andreev current 
peak at h ^ A{T) (for mono-domain case) in the finite tem- 
perature and finite voltage regime, see Fig. [3ja), solid black 
line. However, there are quantitative differences between the 
mono-domain (a = 0) and two domain (a = 7i^ 7i/2) config- 
urations. For a = 0, P{h) shows a minimum at h ^ ^{T). 
It is worth mentioning that around this minimum the cooling 
power of the SIF1F2N system is lower than that of the NIS 
junction (h = 0). By increasing the angle a the minimum is 
less pronounced and shifts to larger values of h> A{T). For 
these values of h and for h = ^ the superconducting conden- 
sate can penetrate both ferromagnetic layers. Thus, the effec- 



FIG. 4. (Color online) Cooling power versus length I2 of the F2 
layer for (a) h = 0.7 A{T) and (b)h=\ .1A{T). We consider different 
orientations of the exchange field vector in the second ferromagnetic 
layer F2 with respect to the one in Fi : a = (solid black line), a = 
71 jl (dashed blue line), a = n (dash-dotted red line), and calculate 
the cooling power at optimal bias; W = 1 x\0~^,T = 0.25Ao and 
/l = ^. P is defined in Fig. 2. 

tive exchange field h acting on the Cooper pairs is a field, aver- 
aged over the length (^/^.^^ The /i(a) is gradually reduced as a 
increase from to TT. As before the cooling power minimum 
is at /i(a) ^ A(r) which in the case of a finite a corresponds 
to larger values of the bare h. The minimum of the cooling 
power (Fig. |2]), corresponds to a maximum of the Andreev 
current [Fig.[3[a)]. The unexpected nonmonotonic behavior of 
the Andreev current at small exchange fields /z ~ A(r) is due 
to the competition between two-particle tunneling processes 
and decoherence mechanisms as quantitatively explained in a 
recent work by the authors. ^^ 

We analyze now the dependence of the cooling power on 
the length of the ferromagnetic bridge F1F2. To do this, we 
fix the thickness of Fi dX l\ = t, and vary I2. Fig. [i] shows 
the P{l2) dependence for two different values of the exchange 
field h/A{T) = 0.7, 1.7 and different magnetic configurations 
a = 0, ;r/2, 71. As expected all curves tend to a finite asymp- 
totic value when h^ ^. This value however depends on a. 

In the case of an exchange field smaller than the supercon- 
ducting gap [h = 0.7A(r), see Fig. |4ja)] the cooling power 
first reduces monotonically to a minimum by increasing I2, 
then enhances to a maximum and finally reduces to the asymp- 
totic value. Such behavior is preserved for all magnetic con- 
figurations and it follows from the nonmonotonic behavior of 
the Andreev current, shown in Fig.|3jb). Decrease of the An- 
dreev current corresponds to the increase of the cooling power 
and vice versa. As shown in Fig.jSjb), at large values of I2 the 
Andreev current increases by decreasing I2, reaches a maxi- 
mum and finally decreases for h^^- The strong suppression 
of the Andreev current for small values of fe is due to the prox- 
imity of the N reservoir at x = /12. On the other hand for larger 
values of I2 the superconducting proximity effect in the ferro- 
magnetic bridge is fully developed and leads to an increase of 
the Andreev current. It is remarkable that the cooling power 
for a = TT is larger than the one at a = for all values of h. 
In this case a lower effective exchange field h leads to larger 
values of the cooling power, due to the shift of the minimum 
of P{h) observed in Fig.|2] 

For an exchange field larger than A(r) [h = 1.7A(r), see 
Fig. |4|b)] the behavior of the cooling power as a function of 
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FIG. 5. (Color online) Cooling power versus bias voltage for h = 
A(r) (a), h=l.lA{T)(b) and h = ^A{T) (c) for different orientations 
of the exchange field vector in the second ferromagnetic layer F2: 
a = (black solid line), a = n/2 (blue dashed line) and a = n (red 
dash-dotted line); W = 1 x 10"^ T = 0.25Ao, h=^ and l2 = 6^.P 
is defined in Fig. 2. 
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FIG. 6. (Color online) Dependence of the cooling power on the tun- 
neling parameter Wfovh = A{T) (a), h=\JA(T) (b) and h = SA{T) 
(c) and for different orientations of the exchange field vector in the 
second ferromagnetic layer F2: a = (black solid line), a = n/l 
(blue dashed line) and a = n (red dash-dotted line). P is calculated 
at optimum bias; T = 0.25Ao, h = ^ and l2 = 6^. We have defined 
P{W,Vopt) = P{W,Vopt)e^Ro/Al. Note the logarithmic scale. 



I2 strongly depends on a. For a mono-domain magnet, a = 0, 
the cooling power increases monotonically by increasing I2 
until it reaches the asymptotic value due to the suppression of 
the Andreev current as in the ballistic case studied in Ref . 1 8 
Similarly, in the antiparallel configuration (a = tt), the cool- 
ing power first increases by increasing I2, however for a larger 
value of I2 reaches a maximum and then decreases. The pres- 
ence of F2 with a magnetization antiparallel to the one of Fi 
leads to a reduced effective exchange field of the F1F2 bridge. 
Thus, the Andreev current contribution is enhanced with re- 
spect to the one in the case /2 = 0. As intuitively expected the 
cooling power (Andreev current) reaches a minimum (maxi- 
mum) when h ^ h = ^, i-e. when the average magnetization 
is minimized. Further increase of I2 > ^ leads to a suppres- 
sion of the Andreev current and therefore to an increase of P 
until the asymptotic values are reached. Fig. |4|b) also shows 
the intermediate case a = 7r/2. 

We now analyze the dependence of the cooling power on 
the bias voltage eV, tunneling parameter W and temperature 
T. In our subsequent analysis we consider three different val- 
ues of the exchange field h = A(r), 1.7A(r), 8A(r), and three 
magnetic configurations a = 0, 7r/2,;r. We set h = ^, short 
enough for the pair correlation to be substantial in the F2 layer 
[for h = A(r), 1.7A(r)] and /2 = 6(^, long enough to ensure 
the asymptotic regime [see Fig.|4|. Figs. [5] [6] and |7] show the 
cooling power as a function of eV,W and T. A common fea- 
ture of these figures is that the range of values of V, W and 



r, for which the cooling power is positive increases by in- 
creasing h. Also the magnitude of the cooling power increases 
with h. This is in agreement with the qualitative predictions 
of Ref. 18 Note that the shape of all curves in Figs. |5j|6] and 
[7] does not depend significantly on the angle a. 

Figs.[5]and|6]show that for low values of eV and W, respec- 
tively, the cooling power depends only weakly on the relative 
magnetization angle a. However, by increasing eV and W the 
difference becomes appreciable, in particular for h^ A{T). 

As shown in Fig. Is] at certain value of eVopt ^ 0.8A(r), 
the cooling power reaches its maximum value Pmax = P^Yopt)- 
The eVopt value is the one used as optimal bias value in the 
figures. For voltages larger than this optimal value, the quasi- 
particle current / and hence the Joule heating power IV in- 
crease drastically leading to a rapid decrease of the cooling 
power. As can be seen from Figs. [5] [6] and [7] the optimal volt- 
age Vopt depends on the temperature T, tunneling parameter 
W and magnetic configuration angle a. For the exchange field 
equal to the superconducting gap the maximal cooling power 
Pmax is largest in the antiparallel configuration, while for larger 
h = 1.7 A(r) the largest value P^ax is in the parallel configu- 
ration, in agreement with Fig.[2j 

Fig. [6] shows that the cooling power has also a maximum 
as a function of W . Increasing W the cooling power first lin- 
early increases as single electron tunneling dominates. For 
larger values of the tunneling parameter, the Andreev current 
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FIG. 7. (Color online) Temperature dependence of the cooling power 
for /i = Aq (a), h = I.VAq (b) and h = SAq (c) and for different 
orientations of the exchange field vector in the second ferromag- 
netic layer F2: a = (black solid line), a = 7i/2 (blue dashed line) 
and a = n (red dash-dotted line). P is calculated at optimum bias; 
W = lx 10-^ li=^ and /2 = 6^. P is defined in Fig. 2. 



heating dominates over the single-particle cooling and leads 
to a rapid decrease of the cooling power, which tends to zero 
at a certain onset point. As the exchange field increases, the 
role of Andreev processes becomes less important, therefore 
the onset shifts towards larger values of W. This means that 
for higher exchange field in the ferromagnetic interlayer one 
may use weaker tunnel barriers for the microcooler fabrica- 
tion, which leads to higher amplitudes of the cooling power 
[see Fig.|6](c)] and more effective electron refrigeration. 

In Fig. [7] we show the temperature dependence of the cool- 
ing power. At r > 0.42Ao ^ 0.757^, where Tc is the criti- 
cal temperature of the superconductor, the cooling power be- 
comes negative for all voltages. This value of the temperature 
holds for a wide range of parameters. ^^ The existence of such 
a maximal temperature is due to the increase of the number 
of thermally excited quasiparticles which produce enhanced 
Joule heat. By lowering the temperature the cooling power 
at optimal bias first increases and reaches a maximum. At 
lower temperatures, the Joule heat due to Andreev processes 
causes the cooling power to decrease. At a certain tempera- 
ture Tfnin, the cooling power tends to zero, which defines the 
lower limiting temperature for the cooling regime. As follows 
from Fig. [7] the temperature Tmin decreases when increasing 
the exchange field; this is because the Andreev current and the 
associated Joule heat are suppressed by the exchange interac- 
tion in the ferromagnet. Finally, one can see from Fig. |7] that 
the minimum cooling temperature in the parallel T^-^ and an- 
tiparallel T^f^ configuration satisfy: T^f^ < T^-^ for h = A(T), 
while T^f^ > T^.^ for h = 1.7A(r). For h = 8A(r) [Fig.JTJc)] 
P{T) is almost independent on a. 

A common feature of Fig. [5] [6] and |7] is that for rather small 



value of the exchange field, h = A{T), the antiparallel con- 
figuration is more favorable for cooling [see (a) panels] . For 
larger exchange field /z = 1 .7A(r) , on the contrary, the parallel 
configuration is favorable for cooling [see (b) panels] . As ex- 
pected, in the case of strong enough ferromagnet [/z = 8A(r)] 
the thickness of Fi layer /i :^ ^^ and the superconducting con- 
densate practically does not penetrate into F2 layer. Thus the 
cooling power is a-independent [see (c) panels of Figs. [5j [6] 
and[7). 



IV. CONCLUSIONS 

We have developed a quantitative theory of charge and heat 
transport in normal metal - superconductor tunnel junctions 
with an intermediate ferromagnetic bilayer. We have assumed 
that the magnetizations of the ferromagnets form an angle a 
and focused our study on the cooling power of such a struc- 
ture. 

In the previous works it has been suggested that the larger 
the exchange field the more efficient the cooling. ^^ In this case 
the enhancement of the cooling is due to the suppression of 
the Andreev processes and therefore suppression of the Joule 
heating, released in the normal metal electrode. However, our 
results have shown that this hypothesis is only valid in the 
case of strong ferromagnets [h:^ A]. For weak ferromagnets 
with an exchange field comparable to the superconducting or- 
der parameter A the cooling power shows a non-monotonic 
dependence on h, with a minimum ath^ A (in mono-domain 
case) that corresponds to a maximum in the Andreev current 
/a. Moreover, around this minimum the cooling power of the 
SIF1F2N structure is even lower than the one of the NIS junc- 
tion. We have also shown that in the two-domain case, a finite 
value of a shifts the minimum of cooling power to larger val- 
ues of h if the thickness of Fi is comparable to the magnetic 
length ^h. In this case, the effective exchange field h acting 
on the Cooper pairs is gradually reduced as a increases from 
to TT. The minimum then is at /i ~ A which corresponds to 
larger values of the bare h. Thus, for exchange fields /z < A the 
antiparallel magnetic configuration (a = ;r) of magnetization 
leads to larger values of the cooling power. Such small ex- 
change fields can be realized in weak ferromagnetic alloy s,!^ 
or in hybrid structures consistin g of fe rromagnetic insulators 
in contact with superconductors .I^^Ell por values of h larger 
than A, the parallel configuration (a = 0) is the one that leads 
to larger values of the cooling power. For values ofh^^A the 
cooling is almost independent of a. 

Finally we have analyzed the dependence of the cooling 
power on the bias voltage, the tunneling parameter and the 
temperature. The optimized values for more efficient cooling 
are shown in Figs. [5] [6] and [7] 
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